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1. Introduction 

Today, the studies of diffusion mechanism and catalytic reactions of adsorbed atoms, formation 
of nanostructures on a metal surface are quite topical in surface physics [1-9]. In the processes of 
adsorption, desorption, and surface diffusion, a metal surface undergoes a reconstruction accom- 
panied by a variation of nonequilibrium properties of both electron and ion subsystems. In this 
case, the electro-diffusive, viscothermal, and electromagnetic properties of the electron subsystem 
change in the field of metal surface ions. Studies of transport processes and their particularities in 
the electron subsystem of semi-bounded metals are of great importance for nanotechnologies and 
catalytic technologies. Generally, the electron subsystem stays in the states far from equilibrium 
during transport processes. This fact significantly complicates the electron subsystem description. 
To discover transport processes, various theoretical approaches are developed for the spatially in- 
homogeneous electrons-atoms systems. Particularly, the time-dependent density functional theory 
(TDDFT) is widely used 13- 1^]. Another theoretical approach is related to the hydrodynamic 



model of surface plasmons for a spatially inhomogeneous electron gas proposed in [23|-|25| with the 
use of the response theory [26] based on the Boltzmann kinetic equation. The quantum statistical 
theory for the description of nonequilibrium processes in the "metal-adsorbate-gas" systems was 
developed in the works using the Zubarev method of nonequilibrium statistical operator 



(NSO) [29, 30]. In particular, a self-consistent description of nonequilibrium processes in the atomic 
and the electron subsystems was presented in j^^j at the kinetic level of the description of electron 
processes. To study the ionic and electron structures of a semi-bounded metal, a generalized ap- 
proach that takes the effect of discreteness of the ion subsystem into account and is based upon 
the model of a semi-bounded "jellium" [si, 32| was proposed in [1, [2^. It is worth noting that the 



effect of discreteness of ionic density on the characteristics of a semi-bounded "jellium" was con- 
sidered in (ssl - fssi by means of constructing a perturbation theory with respect to the electron-ion 
interaction pseudopotential. However, the linear response of the electron subsystem to the lattice 
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potential does not take into account the effects of inhomogeneity of the electron subsystem. The 
approach described in [28,, ,31. . ^2] allows one to model the formation of a surface potential and 
to calculate the partition function for the generalized model in terms of the cumulant averages of 
the "jellium" model. In (isl], the generalized "jellium" model is a basis for the statistical description 
of electro-diffusive processes for the electron subsystem of a semi-bounded metal with the use of 
the NSO method, where the single parameter of the reduced description is the nonequilibrium 
average value of the electron density operator. For such a system, the quasi-equilibrium statistical 
sum was calculated by means of the functional integration method for the case of an electron-ion 
interaction pseudopotential. In principle, it enables us to obtain expressions for the nonequilib- 
rium statistical operator in the Gaussian and in the higher approximations with respect to the 
dynamic electron correlations. In [28], the nonequilibrium statistical operator and the generalized 
transport equation of inhomogeneous diffusion were obtained for weakly nonequilibrium processes 
(linear approximation with respect to the gradient of the electrochemical potential). The same 
approximation is used to find the equation for the "density-density" time correlation function that 
determines the dynamic structural factor of the electron subsystem of a semi-bounded metal and 
to demonstrate the connection of this electro-diffusive model in the linear approximation with 
the TDDFT [10.-13j . The expressions of a nonequilibrium statistical operator in the Gaussian and 
higher approximations with respect to the dynamic electron correlations with the quasi-equilibrium 
partition function calculated by means of the functional integration method were obtained in [s^ . 
It enables us to go beyond the linear approximation with respect to the gradient of electrochemical 
potential. The generalized transport equations for a nonequilibrium average of electron density 
operator for strongly nonequilibrium processes in an electron subsystem of a semi-bounded metal 
were presented for the corresponding approximations of the nonequilibrium statistical operator. 

In this paper we perform a viscoelastic description of the electron subsystem of a semi-bounded 
metal on the basis of the generalized "jellium" model, when the nonequilibrium mean values of 
the electron density operator and the electron momentum density operator are chosen for reduced 
description parameters. These parameters also play an important role in TDDFT [IB, 21|. In the 
second section we find the nonequilibrium statistical operator of a viscoelastic model for an elec- 
tron subsystem of a semi-bounded metal. The general calculation of a quasi-equilibrium partition 
function for a quasi-equilibrium statistical operator of an electron subsystem of a semi-bounded 
metal is obtained using a functional integration method. The found quasi-equilibrium statistical 
operator of the electron subsystem is used for the construction of a nonequilibrium statistical oper- 
ator of the system by means of the Zubarev method (29l . [30| . In the third section we work with the 
Gaussian approximation for a quasi-equilibrium partition function, when the operators of electron 
density and electron momentum density do not correlate as pair correlation functions. Moreover, 
we receive a nonequilibrium statistical operator and the corresponding equations of generalized 
dynamics of the viscoelastic model for nonequilibrium averages of the operators of electron density 
and their momentum density in the Gaussian approximation. In the section 4 we use the following 
higher approximation for a quasi-equilibrium partition function, when static correlations between 
the operators of electron density and electron momentum density occur with the third-order cu- 
mulant averages. In this approximation, the nonequilibrium statistical operator and the non-linear 
dynamics equations of the viscoelastic model of an electron subsystem of a semi-bounded metal is 
obtained. 



2. Nonequilibrium statistical operator for the electron subsystem of a semi- 
bounded metal within the viscoelastic model 

We start with the generalized "jellium" model of a semi-bounded metal that considers the effect 
of the ionic subsystem discreteness. The Hamiltonian of the system could be written as follows: 

H = ^i;a(p)a],(p)a„(p) + ^i^fe(q)pfc(q)p_fc(-q) 

p,a q k 
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q fc q 



§E'-(q|0) + ^ E ^E'^'-(ql^« - e'''^^"'-^"^^ (2.1) 

where — e is the charge of an electron, m stands for the electron mass, r^, z = 1, 2, . . . , iV is the 
electron coordinate, iVion means the number of ions in a metal with relevant charges Ze and 
coordinates Rj-; {—oo< Xj, Yj < +oo, Zj ^ Zq, Zq = const), z = Zq is a dividing surface, V = SL 
is a volume of the system, S means a surface square of a semi-bounded metal, L denotes an area 
of a normal coordinate changing: z G (— L/2,+L/2), S oo, L ^ oo. Ea{p) = + £a is 
the electron energy in the state (p,a). The hachure at sum denotes the absence of terms with 
q = (two-dimension wavevector in a semi-bounded metal surface) due to the electro-neutrality 
condition ZNion = i^fe(q) = Ane"^ / {q^ + k"^) is the Fourier-image of interaction potential between 
electrons, /fc(q) is the Fourier-image of the model pseudopotential of interaction between ions and 
electrons w(ri — Rj) — — |r.^R. + fi^i ~ ^j)j = (^j;^); ^(qN) = 27re^ e~''l^l/g means the 
two-dimension Fourier-image of the Coulomb potential. 



1 



Skiq) = E e-'i""^-'^-^^ (2.2) 

is the structure factor of the ionic subsystem, the Fourier-components of electrons density, the 
"collective" variable is 



p,(q)= E («i|e"=1a2)at^(p)a„Jp-q), (2.3) 

where 



p,Q:i ,a2 



{ai\...\a2) ^ Jdzipl^{z)...ip^^{z), 
ipa{z) and £q, are the eigenfunctions and the eigenvalues of the Schrodinger equation 

d2 



V{r) = V{z) is a surface potential that is a function of only normal electron coordinate. 

The averages of the electrons density operator (/o(r))' and their momentum density (p(r))* could 
be chosen for the main parameters of the reduced description for the study of viscoelastic processes 
within the formulated model. These parameters are connected with the relevant inhomogeneous 
electric and magnetic fields: 

V.(E(r))*^e(p(r))*, 

V X (H(r))' = -i|(E(r))* + --(p(r))*, ^^'^^ 
c at cm 

where c is the speed of light, (. . .)* = Sp[. . . p{t)]^ p{t) means the nonequilibrium statistical operator 
of the generalized "jellium" model, which satisfies the Liouville equation with the Hamiltonian (|2.ip . 
The solution of the Liouville equation for p{t) in the Zubarev method with taking into account a 
projecting technique can be presented in a general form: 

t 

p{t) = p^{t) - j e^(*-*')Tq(t, t') [I - V^')] iLNp^{t')dt\ (2.5) 

— CJO 

where 



T^it,t') = exp+ i J [1 - VS")]iLNdt" 
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denotes the generalized evolution operator with taking into account a projecting technique, Pq{t') 
is the generalized Kawasaki-Gunton projection operator. Its structure depends on the quasi- 
equilibrium statistical operator Pq{t). In our case Vq^t) could be written down as: 



(Pfe(q))* Spp' 



(2.6) 



and it has the following properties 

Kimt) = Pq(i), ^q(i)Pq(i') = Pq(i), rq{t)Vq{t') = Vq{t), (l - Vq{t))Vq{t) = 0. 

We define the quasi-equilibrium statistical operator Pq(t) with the principle of the Gibbs entropy 
maximum when the parameters of the reduced description {p{v)Y and (p(r))* are fixed and the 
normalization condition Sppq(t) = 1 is satisfied. In our case we find 

Pq{t) = exp |^-$(i) - /? jif - ^ [7Ijt(q;t)pfe(q;t) + Afc(q;t) • Pfe(q;t)] | j , (2.7) 

$(t)=lnSpexp ^-/3|H--^^[/Ifc(q;t)pfc(q;t) + Afc(q;t).pfc(q;t)]|j , (2.8) 

where fl). (q; t) = pk (q; t) + e(pk (q; t) stands for the Fourier-component of the electrochemical po- 
tential of electrons, /ife(q;i) is the Fourier-component of the chemical potential, ipk{c{',t) is the 
Fourier-component of the local electrical potential, Ak{q_;t) = Vfe(q;t) — c~^afc(q;t), Vfc(q;t) de- 
notes the Fourier-component of the average velocity of electrons, (q; t) is the Fourier-component 
of the vector potential a(r; t) of electromagnetic field: 



(H(r))* = Vxa(r,i), (E*(r))* 



Id, , 



(2.9) 



where (E*(r))* is a tangential part of electrical field, p^{q;t) is defined with the self-consistent 
condition: 



(Pfe(q))* = W(q))* 



and the thermodynamical relations: 



6^t) 



SS{t) 



S{Pk{ci;t)y 



(Pfe(q))*, 



6S{t) 
5{epk{ci; t))* 



-^Vfe(q;i), 



(2.10) 



(2.11) 



(2.12) 



when (pfc(q;f))* are fixed, where S{t) is the nonequilibrium entropy defined using the Gibbs 
method: 

S{t) = -Sp [In pq{t)]pq{t) 

= m + A{Hr - ;^EE [Mfe(q;OW(q))q + A;.(q;t) • (Pfe(q))*]| 

= lnZ(i)+/3|(i?)*-^^^[7I,(q;i)(pfe(q))'+Afc(q;t)-(pfe(q))*]l, (2.13) 
I q fe J 
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and besides 



(5$(i) 



SSit) 



13 



(v.(q))S 



= (P/c(q))*, 
SS{t) 



afc(q;0- 



is the partition function of the quasi-equihbrium statistical operator: 

Z{t)^Spexp l^-Z^ji/- [7lJq;iK-(q)+Afc(q;t) •pfe(q)] |j 



(2.14) 
(2.15) 

(2.16) 



One should calculate the quasi-equilibrium partition function (j2.16p to define the structure of the 
Kawasaki-Gunton projection operator, to find the nonequilibrium parameters J[f.{ci]t), Afe(q;i), 
and thus, to obtain the nonequilibrium statistical operator p{t). The approach based on the "jel- 
lium" model is proposed in (isl . [s^ in the diffusion description case, where the average value of 
the electron density operator {p{r)y is chosen for the reduced description parameter. We use this 
approach to calculate Z{t) (|2.16p . We apply the functional integration method to present Z{t) as 
follows: 



Z{t) = exp 



ZielAZ{t). 



Here, 



Zjeii^Sp [exp(-/3i/o)T5i(/3)] 



(2.17) 



(2.18) 



is the partition function of the "jellium" model of the electron subsystem that corresponds to the 
equilibrium state, found in [1, [2^ 31. 37|. 



5i(/3) =exp 



I d/3'E'E^^.(q)pfe(q|/?')P-fc(-q|/3') 



q k 



denotes the contribution of electron interaction and Pfe(q|/3') = ^"pfe(q)e 



'P'Ho 



where 



AZit) = — Sp [exp(~/3i?o)T5i(/3)52(/3;0] = (52(/3; i))jcn , 



((• • •))jen = ^ Sp [exp(-/3iJo)T5i(/3)(. . .)] 



(2.19) 



(2.20) 



5*2 (/3;0 = Texp 



p 

^ J d/3'^S(q,fc;Ot^^W(q;/3') 
fc^q 



(2.21) 



and 



C(q, /c; i) = col [Bfc(q;t),Afe(q;t)] 
is a column vector, Bfc(q;t) = NionSkiq)wk{q) -^^(q;^), 

W^i+)(q;/3) = [pfe(q;/3), Pfc(q;/3)] 
is a row vector. One can write down AZ(t) after applying the cumulant representation: 



AZ(t) 



exp 



fcl,...,fe„ qi,...,q„ 

X C'(q„,fc„;t)OT_fci,... _fe„(-qi, . . . ,-q„;/3) 



(2.22) 
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where 



aHfc,,...,fc„(qi, . . . , q„; /3) = i" {TWk, (qi; /?) ■ • ■ T^fc„(q„; /?) 



joll 



(2.23) 



denote matrices of fluctuations of the cumulant averages of electron density and electron momen- 
tum, which are obtained with the nonequilibrium statistical operator of the "jellium" model of the 
electron subsystem d, [28, 31, 32 1. Particularly, the matrix of the second cumulant has the following 
stmctrUrc 

OT^:,fe,(qi,q2;/3) mt^f,,^(qi,q2;/3) 
9?tr,fc,(qi,q2;/3) 971^^,^ (qi, q^; /?) 



9^fci,fc2(qi:q2;/3) 



(2.24) 



where 



2'Jfci,fc,(qi,q2;/3) = (pfei(qi;^)pfe2(q2;/3))jcn - (pfei (qi; ^))jcii(pfe2 (q2; ^))jcii , 



(2.25) 

(2.26) 
(2.27) 



^fcf,fc2('ii'^2;/3) = (pfci(qi;/3)pfc2(q2;^))jon - ((Ofei(qi;/3))jcii(Pfc2(q2;/3))jeii, 
97tfei^fe2(qi,q2;/3) = (Pfci(qi;/3)pfc2(q2;^))jcu - (Pfci(qi;/3))jcii(pfe2(q2;/3))jcii, 

and (pfc; (q^; ^))jeii = 0, since averaging of the momentum density operator proceeds with the 
equilibrium statistical operator. For the same reason, the averages 

(pfci(qi;^)Pfe2(q2;/3))joii (pfci(qi;/3)/Ofc2(q2;/3))joii = o. 



therefore, matrix (|2.24p is a diagonal one: 



9^fei,fe2(qi:q2;/5) 



_/»l^:,fc,(qi,q2;/3) 







mtPP,^(qi,q2;/3) 
Considering the above, one can present AZ{t) in the Gaussian approximation as 



(2.28) 



exp 



Ujz)'^ ^ fi3,,(qi;t)Bfe,(q2;«:,,,(qi,q2;/3) 



ki,k2 qi,q2 



Afc, (qi ; t) ■ ml^j^^ (qi , q2 ; /3) • Ak, (qs ; t) 



(2.29) 



In the next section we find the nonequilibrium statistical operator in the Gaussian approximation, 
where the operators of electron density and electron momentum density do not correlate as a pair 
correlation function. 



3. The Gaussian approximation 



Let us present the relevant statistical operator Pq{t) to find the nonequilibrium statistical 
operator in the Gaussian approximation. Considering (|2.29l) one can write down: 



pf)(t)=exp -5^(0 



where 



5«(i) 



N 



1 V X — ^ 



25* 
q 

1 

2 \SL 



E E [i?fc:(qi;0Sfc2(q2;«^,fc2(qi'q2) 

fci :fe2 qi !q2 

Afc, (qi; • mll^^ (qi, q2) • A^, (qs; t) 

^{^^4zYl + Afc(q; t) ■ Pfc(q)] | 

fc,q ^ 



(3.1) 



(3.2) 
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is the entropy operator. In order to exclude the parameters /i^,(q; t) we use the following thermo- 
dynamical relation: 



from which one can find: 



= (Pfc(q))', 



{PkicDY = ( ^ ) E [5.'(q') -7^.'(q';i)] OT!.'',,,_,(-q',-q), (3.3) 



k' ,q' 



where 5fe(q) = A^io„S'fc(q)wfe(q). 

Denoting [W'']^^^^^(qi,q2) as the inverse function of 971^^ ^.^ (qi, q2): 

E [OT^'']^l.„(q,q")9K^.{'',..'(q",q') = s.^s^,^, (3.4) 

from the Fourier-component of the electrons one can find the electrochemical potential: 



As we can see, the Fourier-component of the electrochemical potential in the Gaussian approxima- 
tion is expressed in terms of the structure factor of the ionic subsystem and the Fourier-component 
of the local pseudopotential of interaction between electrons and ions. Time-dependence is described 
with the average equilibrium value of electrons density, renornialized in terms of the structure fac- 
tor of the ionic subsystem, the pseudopotential Wk (q) and the inverse function [DnPP]-^\{-(i' , -q) 
of the pair irreducible cumulant average of electrons density fluctuations. Similarly we exclude the 
parameters (q; t) from (j3.2p using the thermodynamical relation in the Gaussian approximation: 



SL 

hence, one can find 



<5^Afc(q; t) 



(Pfc(q))* = - ( ^ ) E Afe'(q') • ^ir,i-<i', -q)- (3.6) 



k',q' 

Defining [SCtPP]:^,^.^,, (-q', -q") as the inverse function of (-q', -q"): 

E [mipp]^,i„(q,q")a«^'^,fe'(q",q') = 

k" ,c{" 

for the Fourier-component of A/e(q;t) one can obtain from 

SL 



Afc(q;0 = - ( ^ ) E (Pfc'(q'))* • _,(-q', -q), (3.7) 

k' ,q' 



the time-dependence of which is described with the average value of a density momentum oper- 
ator, renormalized via function [3!Jt'''']l[., _j,(— q',— q), which is the inverse function of the pair 
irreducible cumulant average value of the momentum density fluctuation of electrons. In consider- 
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ation of (|3.5|) and (|3.7p . the entropy operator could be written down as follows: 

= /3^5]'Kq|0)+lnZjen- E E [(p..(qi))*[9Kn:l,,-.,(-qi,-q2)(pfe.(q2))* 
q fei,fe2qiiq2 

+ {Ph(ci,))' ■ [fOTi'-l:;,._(^(-qi. -Cb) ■ {Pb(qj»* 



/c,q 



^ fl- ) EE (p^'(q'))* • [2npp]:l.v-.(-q',-q) • p^(q) 



^fc(q) 



SL 



^(pfc-(q'))'[OT'"']:l.,,_,(-q',-q) 



k' ,q' 



Pfc(q) 



q,q' 



(3.8) 



Then, in order to calculate the nonequilibrium statistical operator (j2.5p in the Gaussian approx- 
imation (|3.ip . (|3.8p one should reveal the structure of the Kawasaki-Gunton projection operator 
(|2.6p . and its effect and the effect of the Liouville operator on pq^-* as it is shown in appendix A. 
Then, considering (|A.3p . (jA.4p for the nonequilibrium statistical operator one can get: 



Pit) 



h'q 



1-r 



W-E /e^<*'-*)rG(M') 

fc,q_oo 

dr ) " (fc, q; t') Lf) • W<^%^ {k,q;t')\ dt' , 



(3.9) 



where 



Jf )(A,q;t') = [l-7'(«)(i')] iipfc(q), I^^H^ci^t') = [l - T^^^^ (i')] iiPfc(q) (3.10) 
are the generalized fluxes, V'^'^^it) is a Mori-like projection operator, which effects as follows: 



p(«)(t)i = EE 

A;,q /i;',q' 



jpfc, (q'; t) [onn :L (-q', -q") Pfc" (q")^ 



<5pfc,(q';t) . [mtpp]:i _,„(-q',-q") • (pfc"(q")^ 



where 



(3.11) 



(...)^ ^ Sp[...p«(i)], 
5pk. (q'; t) = Pk' (q') - (Pfc' (q'))* , -^Pfc, (q'; t) = p^, (q') - (p^, (q'))* ■ 

Using the effect of operator V^'^\t') on iLpfc(q), one can show that Ip'^\k,q;t') = 0. Accord- 
ing to p.9p the nonequilibrium statistical operator in the Gaussian approximation is the func- 
tional of the observable values (pfc(q))*, (pfc(q))* and the generalized fluxes of momentum density 
Ip (fc, q;t'). One can obtain the relevant transport equations for (pfc(q))* and (pfe(q))* with the 
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nonequilibrium statistical operator as follows: 



4(Pfc(q))* + i/cq- • (Pfc(q))* = 0, 
ot m 



(3.12) 



dt 



J2 J e^^*'-*V|,°' (fc, q; k', q'; i, t')W(,«) (fc', q'; t')di' = 0, (3.13) 



where 



l-T 



(3.14) 

is the generalized memory function that describes dissipative processes and takes into account 



(pfc(q))G = - E (Pfc(q)Pfc'(q'))>jp''^(fc',q';0- 

k' ,q' 



(3.15) 



It is the time correlation function of the generalized fluxes of momentum density, averaged with the 
quasi-equihbrium statistical operator p^\t) in the Gaussian approximation. Pfe(q) — iijvPfc(q) = 

ifcq : T)£(q) denotes the tensor operator of a viscous stress of the electron subsystem. This means 
that the generalized memory function p.l4p defines the generalized viscosity coefficient of the 
electron subsystem: 

ip<^p^ {k, q; k', q'; t, t') = -fcq : r/^^ (fc, q; k' , q'; t, t") : fc'q'. (3.16) 



It is notable that the system of equations (j3.12l) , (j3.13p has the same structure as in the case of 
weakly nonequilibrium processes with the only difference in the relevant averages: p^^^ [t) — !• p^{t), 
where 



5pk'{c{,t)[m 



fc',-fc(-q,-q') Jdrp^pkicOPo'" 



1 



<5pfc,(q',t) • [a[n'''']:l„._,(-q,-q') ldTpl-pk{q)Po^ 



where po is the equilibrium statistical operator of the system. 

In the next section we work with the following higher approximation for the quasi-equilibrium 
partition function, when the static correlations between the operators of electron density and 
electron momentum density occur with the third-order cumulant averages. 

4. The quasi-equilibrium partition function in the approximation of the third- 
order cumulant averages 

According to (j2.17p . (j2.22p the relevant statistical operator Pq{t) in this approximation is 



p,«+i(t)=exp f-5(«+i)(i) 



(4.1) 
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<i,(G+l)(i) = 



1 / /3 



2 V S-i 



E E [Sfe:(qi;«^,fc,(qi,q2)i?fc.(q2;i) 



*:i,fc2 qi.q2 



Afe, (qi; • OTPPfc^ (qi, qa) • A^, (qa; 



^(^) E E [i?/c,(qi;t)Sfe.(q2;t)i3fc3(q3;«::'fc„fc3(qi,q2, 



qsj 



kiMM qi:q2:q3 
rppp 



3Afc, (qi ; • mtPP^^ (qi , qa , qa ) • A^, (qa ; t)Bk, (qa ; 



(4.3) 



In order to obtain an explicit form of the nonequilibrium statistical operator and the transport 
equations within the approximation (|4.1[) one should exclude parameters /Ij, (q; t) and A^ (q; t) 
from (|4.2p . In the same way as in the Gaussian approximation case we use the thermodynamical 
relations: 



i5A7ife(q;i) 



(4.4) 



,5$(g+i)(t) 
5^Afc(q;t) 



(Pfe(q))*- 



(4.5) 



Considering the structure of '^'^'^^^\t) in S'('^+^^(t) from (|4.4p one can get equations to define the 
electrochemical potential of the electron subsystem of a semi-bounded metal within the generalized 
"jellium" model: 



(pfc(q))* 



SL 



E Sfe,(qi;t)OT^;',fc(qi,q) 



fei.qi 
i / /? 



M^j E E [sfe.(qi;«:a-2,.(qi'q2,q)Sfc2(q2;0 



*:i.'£2 qi,q2 
>ppp 



3Afe, (qi; t) • SJIPP^^ Jqi, q2, q) • Afe, (qa; t) 



(4.6) 



and to define Afc(q; t) from (|4.5p one can find: 



Pfc(q))' - E Afc.(qi;i)-2?i^,^(qi,q) 



fei.qi 
2 



i(^) E E A..(qi;0-9?t^rfc2,/c(qi:q2,q)s.2(q2;t). (4.7) 



ki,k2 qi,q2 



The detailed calculations are presented in the appendix B. Considering (IB.12P the nonequilibrium 
statistical operator of the electron subsystem of a semi-bounded metal within the generalized 
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"jellium" model, according to (|2.5|) we present as follows: 



X < / dr 





dr 



(4.8) 



where 



Pfc(q) = iLArPfc(q), Pfc(q) = iLArpfc(q). 



(4.9) 



Then, taking into account the above found we obtain the transport equations for the reduced 
description parameters (pfc(q))*, (pfe(q))*: 

t 

|(p,(q))' = (p.(q))*G+i)-E /e^^*'^*)<i>^r'^(^'q'fc''q';^'0<''^'^(fc'q;*')di' 



^ J e^(*'-*)$(G+i) (fc, q, fc', q'; t, t') ■ W(';+i) (A, q; t')At' , 



(4.10) 



^(P.(q))* 



(pfc(q))(G+i) - E /e^^*'"'^<p+'^(fc,q,fc',q';i,i')w^ip''+'^(fc,q;Odi' 

t 

^ /'e^(''-')$(';+i)(fc,q,fc',q';i,t') : W(<;+i)(fc,q;i')di', (4.11) 



where 



((...))jG+i)=Sp(...)pf+i)(i), 



4^+'^ (fc, q, fc', q'; t, t'), 4°+'^ (fc, q, fc', q'; t, t'), ^^^^'^ {k, q, fc', q'; i'), *pp^'^ {k, q, fc', q'; *') 
are the generalized transport kernels that describe diffusive, visco-diffusive and viscous processes 
of the electron subsystem of a semi-bounded metal within the generalized "jellium" model. Partic- 
ularly, ^pp^^"* (fc, q, fc', q'; t') and 'I'pp^^'' {k, q, fc', q'; t') have the following structure: 

$(^+l)(fc,q,fc',q';t,i') = (pfe(q)rf (i, t') fl - (i') 



dr 



pf+i)(Ol p.^(q')U°+^'(0 



(G+1) 



fcq-7:>(°+l)(fc,q, fc',q';i,t') -fc'q', 



(4.12) 



<i>(<;+i)(fc,q,fc',qV,0 = (p.(q)Tf+l)(t,0fl-7'f+i)(0 



dr 



pf+i)(0 p,,(q') 



-r \ t' 

(G-t-l) 



fcq:77(G+i)(fc,q,fc',q';i,0 :fc'q', 



(4.13) 
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1?^'^+^) (fc, q, fc', q'; t, t'), 77'*^+^^ (fc, q, fc', q'; t') denote the coefficients of nonlinear diffusion and 
nonlinear viscosity of the electron subsystem of a semi-bounded metal within the "jellium" model. 
The generalized transport equations (|4.10p . (|4.1ip and the nonequilibrium statistical operator are 
strongly nonlinear ones in comparison with the transport equations (j3.12p , p.l3p , which correspond 
to the Gaussian approximation for pq (t). 



5. Conclusion 



Viscoelastic processes in the electron subsystem of a semi-bounded metal are described on the 
basis of the generalized "jellium" model with the use of the NSO method, where the parameters of 
the reduced description are the nonequilibrium average values of electron density and electron mo- 
mentum. Applying the functional integration technique, we have calculated the quasi-equilibrium 
partition function for such a system in the case of the model pseudopotential of electron-ion inter- 
action in a metal in the Gaussian approximation and in the following higher approximation, where 
the static correlations between the operators of electron density and electron momentum density 
are taken into account with the third-order cumulant averages. 

We have also obtained the expressions for a nonequilibrium statistical operator, which enables 
us to go beyond the linear approximation with respect to the gradients of electrochemical potential 
and average electron density. In the respective approximations for a nonequilibrium statistical op- 
erator we have derived the generalized transport equations for a nonequilibrium average values of 
the electron density operator and the electron momentum density operator that can be applied to 
the description of strongly nonequilibrium processes for an electron subsystem of a semi-bounded 
metal. The generalized transport coefficients (that related, for example, with the generalized vis- 
cosity coefficient of the electron subsystem of a semi-bounded metal) which are contained in the 
corresponding transport equations are calculated using the quasi-equilibrium statistical operator 
in the respective approximations: the Gaussian one (13.11) and the following higher approxima- 
tion (|4.ip . An important point in such an approach is that the time correlation functions and 
the generalized transport coefficients are calculated with the quasi-equilibrium statistical opera- 
tor in the corresponding approximation and represent the functionals of the observable quantities 
{Pk{(l)y, (Pfc(q)}* of a certain order. Of special interest in this approach are the investigations of 
the dynamic structure factor for the nonequilibrium electron subsystem of a semi-bounded metal. 

A. Kawasaki-Gunton projection operator in the Gaussian 
approximation for p^^\t) 

Considering ^J^, ([SH) and 

[$Hn:l.'.-fc(-q',-q)pfHi), (a.i) 



s (pfc(q))* 

<5(Pfe(q))* 



E 

k' ,q' 

E 

k' ,q' 



Pfe,(q';T;t)-(p,,(q'))* 



where 



Pfc,(q';r;i)= / dr p'^'^^t) p,,(q') pf\t) 



Pfc,(q';r;t)= / dr pf\t) Pfc,(q') pf\t) 
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we find the Kawasaki-Gunton projection operator: 



X [aJin:^,-fe(-q', -q) {pM)f Pf\t) Sp(p') 

+E E (p'^"(qV'0 - (Pfe'(q'))*) • _,(-q',-q) • pf \t)Sv{p') 



-EE [9«'"']:l.',-fe(-q',-q)Sp(Pfe(q)p')pfni) 

A;,q fc',q' 

-EE (pfe'(q';T,t) - (p..'(q'))*) • [9JtPP]:i, _,(-q', -q) • Sp(pfc(q)p') i W(A.3) 
Taking into account: 

1 

■^LNpf\t) = 5^P^(,G)(A,q;0 /drLf)(d'p,(q)UG)(t) 



fc,q 



^ W(<;)(fc,q;i) y dr [pf )(t)J • pfe(q) [pf )(0 



(A.4) 



where 



^^^"^^'iiz) E (Pfe'(q'))*[9:«ni,-.(-q',-q) 



A;' ,q' 



W(';)(fc,q;t) = ^(p,,(q'))*-[2nPP]:^__,(-q',-q). 



(A.5) 
(A.6) 



/c' ,q' 



B. Approximate determination of /^^(q; i) in p^^\t) 

Then, we apply the approximation [36] that linearize the equations (j4.6p . (|4.7p with values 
/ifc(q;i) p.6p and Afc(q;t) (|3.8p in the Gaussian approximation. So in (|4.7p one can obtain: 



(pfc(q))* 



SL 



J2 Afe,(qi;i)-Gg,(qi,q;i), 



(B.l) 



fci,qi 



Gr^yqi,q;t) = 2nPP,(qi,q) 



i 5] (q',q2)5mPP^,,,(qi,q2,q). (B.2) 

k' ,k2 q%q2 



Defining the inverse function [G^'^^j ^ (qi, q; t) of G^'^''^(qi, q; t): 



E [g^^^ 

fc",q" 



(q', q"; OGi"^,fc,(q", qi; = -^fc^fc, V.qi 



we find from (jB.ip the following: 



Afci(qi;0 = - 



A;',q'' 



(B.3) 
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Considering (jB.3|) with the Hnearization (|4.6p , we obtain the equation for defining /ij. (q; t) : 



(P/c(q)) 



^(5fe,(q')-7I,-(q';t))9Jl^f,(q',q) 



SL 

+ ^JZ ^ ^ (pfe'(qO)*[»t'''']fe^feJq^ql)(^fe.(q2)-Mfc,(q2;^))OT^:a,,,(ql,^ 

ki,k2,k' qi,q2,q' 

I E E (Pfc'(q'))* • KPP,^.,,(q',q",q) • (p."(q"))* , 



k' ,q' 



k' ,k" q',q" 



(B.4) 



where 



K^F^,.(q',q",q)= E E [^^'^ 

fcl,qi fc2,q2 



_ (q',qi;i«f.^,,,,.(qi,q2,q)[G(3)l" (q2,q";i) 

k' .ki ' \- A k2 -k" 



(B.5) 

is the function that talces into account the complicated renormahzation of the cumulant average 
"momentum-momentum-density". Then, selecting terms from /ife(q;i), we rewrite (jB.4p : 



(pfc(q))* 



^GW(q) - ^ E M.'(q';i)Gi'!.(q',q;0 



k' 



3i 
2! 



E E (Pfc'(q'))*-KPP,^^,(q',q",q)-(Pfc"(q"))* 



(B.6) 



k' ,k" q',q" 



where 



Gi'^(q) = E [^^-(q')»t^f'.,(q',q) 



fc' ,q' 



+ iE E ('^'^■'(q'))*[2^'"]^'!fc.(q'>qi)^fe.(q2)sH^::'fe2,fe(qi'q2'q)]' (b-t) 

fci,qi fe2,q2 



Gi'lfc(q'>q;0 = wt^^(q',q) 



+ E E (^^'(q'))*[»^'"']fc'k(q''qi)»^rfe2,fe(qi'q2.q) 

fel.qi fc2,q2 



(B.8) 



Defining [G*^^'] ^ (q2, q; i), the inverse function of G^i^,\{q' , q; t): 

E f^^'l'' (q',q";0Gi''\,(q".qi;i) = '5fc',fciV, 

^ — ' L J k' .k" ' 



',qi ' 



fc",q" 



one can get from (jB.6l) : 

Aifc(q;0 = - 



k' ,q' 



{Pk'{c^))'-§Y^G^}{c^) 



(q',q;i) 



3i (£_ 
2! I S-L 



E E (p^^'(q'))*-Kr-^.'",.(q',q",q;0-(p."(q"))* 

k' ^k" q',q" 



where 



Kr,pr„,,(q',q",q;0= E K^'T,^. (q'' q"' qi^ ^) G^'^ ^ , (qi'q;^)- 

Kl,qi 



(B.9) 



(B.IO) 
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Now considering (IB.9|) . (|B.3p we can write down the entropy operator S''^'^+^^(i) (14. 2p as follows: 



EE 

3i 



(pfe'(q')>*-^G«(q') 



(q',q;i)pfc(q) 



I E E (Pfc'(q'))* • Kr,^r",.(q',q",q;i) • (p."(q"))*Pfc(q) - f^^+^^W, (b.h) 



k' .k" q' ,q'' 



where according to (|4.3p and (jB.Sp . (jB.9p is a function of a higher order with respect 

to the parameters of the reduced description (pfc(q))*, (pfc(q))*, including the sixth with respect 
to (pfc(q)}* and its structure is important in calculating the Kawasaki-Gunton projection opera- 

tor (12. 6p in terms of —2 — . — -3 — — — -— . Considering the above, we get 



1 



5:ydr[pf+i)(t)] p.(q) 



{t) " .W(';+i)(fe,q;0, (B.12) 



where 



W-(«+i)(fc,q;0 = ^ 



W(G+i)(fc,q;i) 



3i 



(pfc'(q'))*-^GW(q') 



fc'.fc 



(q',q;i) 



:^ E E (p.'(q'))'-Krr^,fc(q',q",q;i)-(Pfc"(q"))*, (B.13) 



2! 

k' ,k" q' ,q'' 

E(p^'(q'))' 

fc' ,q' 



G(3) 



k',k 



(q',q;i)- 



(B.14) 
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B'si3KO-ejiacTi/iHHMM oni/ic HepiBHOBa>KHoY e/ieKTpoHHoY niflcucreMi/i 

HaniBo6Me>KeHoro Mera/iy 

B yaara/ibHeHiw MOfle/ii ">Ke/ie" 

n. n. KocTpo6iiP, B. M. MapKOBi/iJil, A. I. Baci/i/ieHKcP, M. B. ToKapnyhlil^l 

HaL4ioHa/ibHw£i yHiBepcwTer "/IbBiscbKa ncnijexHiKa", Byji. C. BaHflepw, 79013 /IbsiB, VKpaTHa 
iHCTMxyT 4>i3WKW KOHfleHCOBa HMX CHCTeM HAH yKpaiHi/i, Byji. I. CBeHL4iL4bKoro, 1, 79011 JlbBiB, VKpaiHa 

3anponoHOBaHO a'aBKO-e/iacTUHHi/ii/i oni/ic ejieKTpoHHoT niflCHcreMU HaniBo6Me>KeHoro Mexajiy Ha ocHosi ysa- 
rajibHeHoT MOfle/ii ">Ke/ie" is BacxocyBaHHJiM Mexofly HepiBHOBa>KHoro cxaxMCXMHHoro onepaxopa 3y5apeBa. 
OxpuMaHO HepiBHOBa>KHMM cxaxMCXHHHHH onepaxop xa BiflnosiflHi ysarajibHeHi piBH^HHa nepeHocy fljia He- 
piBHOBa>KHMX cepeflHix BHaneHb onepaxopiB rycxi/iH Hi/ic/ia ejieKxpoHis xa Tx iMny/ibcy y raycoBOMy xa Bi/iinoMy 
3a HUM Ha5/iM>KeHi, mo BiflnoBiflae KyMy/i5iHXHi/iM cepeflHiM xpexboro nopaflKy npn pospaxyHKy KBasipiBHO- 
Ba>KHoT cxaxwcxuHHoT cyMM MexoflOM 4>yHKL4ioHajibHoro iHxerpyBaHHa. 

KyiKJHOBi c/iOBa: ysarajibHena MOftenb ">Ke/ie", HepiBHOsa/KHnit CTaTi/icTi/iHHi/ii^ oneparop 3y6apeBa, 
HaniBo6Me>KeHiAi^ Meraji, pisHfiHHfi nepenocy, KBa3ipiBHOBa>KHa CTaTi/icTiAHHa cyma 
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